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The cumulative (anti)neutrino production from all core-collapse supernovae within our cosmic horizon gives 
rise to the diffuse supernova neutrino background (DSNB), which is on the verge of detectability. The observed 
flux depends on supernova physics, but also on the cosmic history of supernova explosions; currently, the 
cosmic supernova rate introduces a substantial (±40%) uncertainty, largely through its absolute normalization. 
However, a new class of wide-field, repeated-scan (synoptic) optical sky surveys is coming online, and will map 
the sky in the time domain with unprecedented depth, completeness, and dynamic range. We show that these 
surveys will obtain the cosmic supernova rate by direct counting, in an unbiased way and with high statistics, 
and thus will allow for precise predictions of the DSNB. Upcoming sky surveys will substantially reduce the 
uncertainties in the DSNB source history to an anticipated ±5% that is dominated by systematics, so that the 
observed high-energy flux thus will test supernova neutrino physics. The portion of the universe {z < 1) 
accessible to upcoming sky surveys includes the progenitors of a large fraction (~ 87%) of the expected 10 - 
26 MeV DSNB event rate. We show that precision determination of the (optically detected) cosmic supernova 
history will also make the DSNB into a strong probe of an extra flux of neutrinos from optically invisible 
supernovae, which may be unseen either due to unexpected large dust obscuration in host galaxies, or because 
some core-collapse events proceed directly to black hole formation and fail to give an optical outburst. 

PACS numbers: 98.80.-k, 97.60.Bw, 95.85.Ry, 98.70.Vc 



I. INTRODUCTION 

Core-collapse supernovae are the spectacular outcome of 
the violent deaths of massive stars. These events, which in- 
clude Type II, Type lb, and Type Ic supernovae, are in a real 
sense "neutrino bombs" in which the production and emission 
of neutrinos dominates the dynamics and energetics. This ba- 
sic picture now rests on firm observational footing in light of 
the detection of neutrinos from SN 1987A iQl, U. Thus all 
massive star deaths - certainly those that yield optical explo- 
sions, and even "invisible" events that do not - are powerful 
neutrino sources. Yet only the very closest events can be in- 
dividually detected by neutrino observatories, leading to burst 
rates so small that no new events have been seen in more than 
two decades. 

All core-collapse events within the observable universe 
emit neutrinos whose ensemble constitutes the diffuse super- 
nova neutrino background (DSNB) ^ [3-10]. Core-collapse 
supernovae produce all three active neutrino species (and their 
antineutrinos), all in roughly equal numbers. However, for 
the foreseeable future only the v^. flux can be detected above 
backgrounds present on Earth. Specifically, the DSNB dom- 



Type la supernovae do not have substantial neutrino emission > 10 MeV, 
but an intriguing alternative fate of accreting wliite dwarfs is the accretion- 
induced collapse (AIC) to a neutron star. Ref . 1 11 ] suggests the AIC events 
can also produce neutrino emission similar' to core-collapse events. If so, 
AIC events would contribute to the DSNB and to optically visible out- 
bursts. However, these AIC events have not yet been observationally con- 
firmed and the expected AIC rate is much lower than that of core-collapse 
events. Therefore AIC neutrinos should not greatly change our results. 



inates the (anti)neutrino flux at Earth in the ^ 10 — 26 MeV 
energy range, and has long been a tantalizing signal that has 
become a topic of intense interest (e.g., lIs UlOll ). Until now 
no DSNB signal has been detected, which set an upper bound 
on the DSNB flux. Super-Kamiokande (Super-K) set the up- 
per limit to be 1.2 cm^^ s^^ above 19.3 MeV of the neutrino 
energy fl5]. However, this limit is akeady close to theoreti- 
cal prediction and thus Super-K is expecting to detect the first 
DSNB signal within the next several years. 

Recently, Ref. [10] considered a variety of complementary 
indicators of the cosmic supernova rate, and concluded that 
the DSNB is no more than a factor --2-4 below the 2003 
Super-K limit [12.1 . Moreover these authors point out that if 
Super-K is enhanced with gadolinium to tag detector back- 
ground events (v^, the resulting enhanced sensitivity at 10 - 
18 MeV should lead to a firm DSNB detection. 

In light of the impending DSNB detection it is imperative 
to quantify the uncertainties in the prediction and to reduce 
these as much as possible. The predicted flux depends cru- 
cially on: (a) supernova neutrino physics, via the emission 
per supernova; and (b) the cosmic history of core-collapse su- 
pernovae, via the cosmic supernova rate (hereafter, CSNR). 
Our emphasis in the present paper is on the CSNR, which has 
begun to be measured in a qualitatively new way by "synop- 
tic" surveys. These new campaigns repeatedly scan the sky 
with a certain fields of view and high sensitivity. Pioneer- 
ing synoptic surveys are already in hand and have shown the 
power of this technique. To date, these surveys have reported 
the detection of several hundreds of supernovae in total, in- 
cluding both Type la and core-collapse events [14-21]. Fu- 
ture surveys, such as DES, Pan-STARRS, and LSST, should 
find > 10"' CCSNe yr^^, eventually with detection rates of 
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> 10^ CCSNe yr~^ based on their depths and large fields of 
view 12211 . Current predictions show that these data will pro- 
vide an absolute measurement of the CSNR to high statistical 
precision out to z ~ 1 ll22i l231. Note that observations seem 
to suggest that Type Iln supernovae are intrinsically the most 
luminous core-collapse type ll24ll . and therefore would con- 
tribute to most of the detections at z > 0.5; but as we discuss 
below, the nature of the bright end of the supernova luminos- 
ity function remains uncertain and other rare but bright super- 
nova types ifTsl I25I - I27I1 might also be important at these large 
redshifts. 

It is important to appreciate that the most crucial input 
from future synoptic surveys will be the normalization of the 
CSNR. The shape of the CSNR follows from that of the star- 
formation rate due to the very short lifetimes of all massive 
star progenitors, and the cosmic star-formation redshift his- 
tory is already relatively well-known out to z ^ 1. However, 
the CSNR normalization is only known to within ^ 40%. 
This will be greatly improved by future synoptic surveys, 
which should measure the CSNR to extremely high precision 
at z ^ 0.3, and therefore dramatically reduce the uncertainties 
in the CSNR (and hence the DSNB) normalization. 

Because our focus is on the interplay between synoptic sur- 
veys and neutrino observations, we wish to carefully distin- 
guish different outcomes for massive stars and their result- 
ing optical and neutrino emission. All collapse events pro- 
duce neutrinos; however, simulations have shown that both 
the amount and energies of the supernova neutrinos varies 
with the mass range of the progenitor stars and how they end 
their lives fG, 10, 28- 32]. Unfortunately, there exists great un- 
certainty about the fate of massive stars, and the as-yet unre- 
solved physics of the baryonic explosion mechanism rnay well 
pla y an important role in determining the outcomes Jl^, l36l - 
Qy]. Recent work suggests that stars below some characteris- 
tic mass (estimated at ^ 25Mq) do explode, producing opti- 
cal supernovae and leaving behind neutron stars; on the other 
hand, stars above some mass scale (estimated at ^ 40A/o) are 
expected to collapse directly into massive black holes without 
optical signals 0281 l3ll l33(l . It is possible that between these 
regimes, a mass range exists (e.g., 25 — 4OAf0) that would 
be a gray area where core-collapses form black holes from 
fallback while still being able to display some (perhaps dim) 
optical signals. 

In the following sections, we will refer to those massive 
stars that first undergo regular core collapse and bounce as 
"core-collapse" events, whether they ultimately leave behind 
neutron stars or black holes formed from fallback. Those 
massive stars that collapse directly to black holes we will re- 
fer to as "direct-collapse" events. Events that also produce 
substantial electromagnetic outbursts we refer to as "visible"; 
those that do not are "invisible." For simplicity, but also fol- 
lowing current thinking, we take visible events to be core- 
collapse events that produce neutron stars and conventional 
(i.e., 1987A-like) neutrino signals. We take invisible events to 
be direct-collapse events, which have a higher-energy neutrino 
signal [31]. "Failed" supernovae should be invisible from our 
viewpoint, though some may have weak electromagnetic sig- 
nals that we henceforth ignore [34]. 

The focus of this paper is to quantify how the CSNR de- 



termination by future synoptic sky surveys will improve the 
DSNB prediction, and to point out some of the science pay- 
off of this improvement. After summarizing the DSNB cal- 
culations {^^, we present our forecasts for the CSNR mea- 
surements by synoptic sky surveys ( ^TUi . Using these, we 
show the impact on the DSNB ( ijlVI) . In particular, we dis- 
cuss present constraints on invisible events, and strategies for 
DSNB data to probe the fraction of massive star deaths that 
are invisible (^Vji. We then switch to a extremely conserva- 
tive viewpoint and discuss the robust lower limit on the DSNB 
( ijVIl) . Conclusions are summarized in Will 



II. DSNB FORMALISM AND PHYSICS INPUTS 

The neutrino signal from the ensemble of cosmic collapse 
events is conceptually simple, and is given by the line-of-sight 
integral of sources out to the cosmic horizon (more precisely, 
to the redshift where star formation begins; in practice, the 
result does not change once redshifts of a few are reached). 
The well-known result is 
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where (e) is the neutrino intensity (flux per solid angle) of 
cosmic supernova neutrinos with observed energy e. Because 
Earth is transparent to neutrinos, detectors see a total (angle- 
integrated) flux (e) from the full sky. Note that neutrinos 
and their energies are measured individually, so the intensity 
and fluxes measures particle number, not the energy carried 
by the particles. Two source terms, 7?,tot and Nv[{\ + z)e], 
appear in Eq. [T] 7?.tot is the cosmic rate of collapse events, 
i.e., the number of collapse events per comoving volume per 
unit time in the rest frame. Each source, i.e., each collapse 
event, has a neutrino energy spectrum A'^i^(eomit) in its emis- 
sion frame with rest-frame energy eomit = (1 + z)e; the factor 
(1 + z) accounts for the redshifting of energy into the ob- 
server's frame. Because we allow for different neutrino en- 
ergy spectra for core-collapse (CC) and direct-collapse (DC) 
events, Nv[{l + z)e] can be expressed as 



N,[{l+z)e\ = /cciV. [(l+2)e]+/DC (2) 

where /dc = '^to?/'^tot and fee = 1-/dc are flie fractions 
for direct-collapse and core-collapse events, respectively; we 
assume these to be constants independent of time and thus 
redshift. Because these fractions are very uncertain, below 
we will consider a range of possible values. Finally, for the 
standard ACDM cosmology the time interval per unit redshift 
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Equation [T] shows that three inputs control the DSNB: (i) 
cosmology, via the cosmic line integral and parameters; (ii) 
supernova neutrino physics, via the source spectrum, (iii) as- 
trophysics, via the CSNR. Of these, the cosmological inputs 
entering via Eq. [3] are very well understood and their error 
budget is negligible. We adopt the standard ACDM model, 
with parameters from the 5-year WMAP data: flm = 0.274, 
Ha ^ 0.726, and Hq = 70.5 km Mpc'^ S. Within 
this fixed cosmology, DSNB predictions require knowledge 
of the source spectra and CSNR. The purpose of this paper is 
to forecast the effects of future improvements on the source 
rate, but to illustrate these we must adopt source spectra. 

Core-collapse neutrino spectra are in principle calculable 
from detailed supernova simulations, e.g., ll29l [36l - l39ll . In 
practice, it remains quite difficult to simulate supernova neu- 
trino emission accurately within realistic explosion models (if 
they explode at all!) and certainly it remains computationally 
prohibitive to perform such ab initio simulations over wide 
ranges of supernova progenitors. Consequently, in DSNB pre- 
dictions different groups have taken different approaches in 
estimating neutrino energy source spectra. Here, we adopt the 
treatment in the recent DSNB forecasts of Ref. [IC^]. These 
authors approximated the neutrino energy spectra as Fermi- 
Dirac distributions with zero chemical potential: 
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where £^ is the total energy carried in the electron antineu- 
trino flavor and is the effective electron antineutrino tem- 
perature. Neutrino flavor change effects are absorbed into 
the choices of £^ and T^. Following Ref. lilOJ . we assume 
the total energy is equally partitioned between each neutrino 
flavor for both core-collapse and direct-collapse events, i.e. 
El, = f i/,tot / 6 for individual neutrino flavor, where £i,,tot is 
the total (all-species) energy output. The variation in neu- 
trino emission from different core-collapse progenitor stars is 
in general expected to be small because neutrinos come from 
newly-formed neutron stars. We adopt £ij,tot = 3 x 10^"^ erg 
per core-collapse event. Ref. |10] finds that the average tem- 
perature after neutrino mixing is constrained to lie in the range 
T,, 4 - 8 MeV. We choose = 4 McV as our benchmai-k 
temperature, which is close to the empirically-derived spec- 
trum ofSN 1987A ||40ll. 

For the direct-collapse events, hydrodynamic simulations 
show that the neutrino spectra are sensitive to the progeni- 
tor masses and nuclear equation of states, with models giv- 
ing total neutrino energy outputs ranging from 1.31 x 10^'^ to 
5.15 X 10^"^ erg and different neutrino average energies rang- 
ing from e^^s = 18.6 to 23.6 MeV |l29l43I|]. We choose the 
model with higher energy so it will create a greater difference 
for comparison. That is, we take £,y,tot — 5.2 x lO'^'^ erg, and 
Tp^ = e^^s/3-15 = 7.5MeV. 

In what follows, we first take all supernovae to be core- 
collapse events (thus visible) as the fiducial case, and then we 
will examine the impact of the direct collapse (invisible) su- 
pernova scenario. Since the emission from the direct-collapse 
events is taken to be larger, this will increase the DSNB de- 
tection rates. Cosmic supernova neutrinos will be detected 



mainly via inverse beta decay iJ^ + p —i' n + interactions 
with protons in a liquid water or scintillator detector. This 
reaction is endoergic with the threshold energy of 1.8 MeV. 
To a good approximation, the nucleon remains at rest, so that 
eg+ ~ e — A, where is the positron total energy, e is the 



i^e energy, and A = m„ 



1.295 MeV. The expected 



differential event rate, per unit time and energy, is 
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The well-known inverse beta decay cross section a^pie) 
ll4ll I42I1 . taken here at lowest order, and which increases 
with energy roughly as Thus the event rates give larger 
weight to the high-energy neutrino flux, which, as we will see 
is the regime best probed by supernova surveys. The total 
event rate in a detector sensitive to neutrino energies e is thus 
R ^ /^'■"*'' dR/de de. The factor Np in Eq.|5] gives the num- 
ber of free protons (those in hydrogen atoms) in the detector; 
in our calculations, we use the value corresponding to 22.5 
kton of pure water for Super-K. 

The upper panel of Fig. [T] shows the neutrino event rate - 
the integrand of Eq. (1) with T,^ = A MeV - with respect to 
redshift at certain fixed observed energies. Because of red- 
shift, neutrinos with low observed energies are more likely 
to come from high redshift supernovae, while neutrinos with 
high observed energies are more likely to come from low red- 
shift supernovae. 

A measurement of DSNB neutrinos and their energy spec- 
trum will thus provide unique new insights into the physics 
of massive-star death. But for the DSNB to usefully probe the 
neutrino emission from supernova interiors, the cosmic source 
rates must be known. It is to this that we now turn. 



III. DSNB ASTROPHYSICS INPUT 

The CSNR not only controls the DSNB flux, but also is of 
great intrinsic interest, and has a direct impact on numerous 
problems in cosmology and particle astrophysics. The stellar 
progenitors of both core-collapse and direct-collapse events 
are very short-lived; consequently the CSNR is closely related 
to the cosmic star-formation rate, which has been intensively 
studied for the past decade Il43i1 . From the present epoch back 
to z ~ 1, the cosmic star-formation rate increases by an order 
of magnitude. At higher redshifts, z > 1, the cosmic star- 
formation rate becomes less certain, but the z < 1 regime 
is responsible for a large fraction of the observable DSNB 
signal. On the other hand, while the shape of the cosmic 
star-formation rate is relatively secure, the absolute normal- 
ization remains harder to pin down. Recent estimates using 
multiwavelength proxies for the star-formation rate indicate a 
±20% uncertainty at z = and a larger uncertainty at higher 
redshift, producing an average of ±40% uncertainty on the 
DSNB detection rate 1 10]. For the direct supernova rate data 
reported in Ref. llOj] , here we adopt a ±40% uncertainty at 
z = 0, double that on the star-formation rate itself (this should 
not be confused with the 40% above). 

Fortunately, a new generation of powerful sky surveys are 
poised to offer a new, high-statistics measure of the CSNR. 
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FIG. 1: DSNB and synoptic survey redshift distributions. Upper- 
Panel: The integrand of Eq.[T]as a function of redshift for different 
choices of observed neutrino energies; this shows the redshift distri- 
bution of sources that contribute to the DSNB signal at these ener- 
gies. Here we assume all the supemovae are core-collapse events, 
as defined in Bottom Panel: The blue curve is the supernova de- 
tection rate by LSST in r-band as a function of redshift, with survey 
depth mf^^ = 23""'"5 and sky coverage of 6.1 sr (20,000 deg^). The 
black curve is a more conservative estimation of the LSST supernova 
detection rate by excluding Type Iln supemovae, which seem to be 
the most luminous based on the small sample of current data. The 
red curve is the fiducial supernova rate for comparison, which is the 
full-sky supernova rate without considering dust extinction or survey 
depth. The curves have bin size Az = 0.05, and the band thickness 
(which are in most cases thinner than the curve width) represent the 
statistical uncertainty 1 / vTV. 



These surveys have wide fields of view and large collecting 
areas, in order to produce deep scans of large portions of 
the sky. These synoptic surveys are designed to repeatedly 
scan a large portion of the sky every few nights with limit- 
ing single-exposure magnitudes of ^ 21'°''^ to ^ 24'"''^, and 
possibly deeper in several passbands. Relatively more mod- 
est prototype synoptic surveys have akeady been completed, 
e.g., SDSS-II ill and SNLS 1 16], or are underway, e.g., the 
Pan-STARRS 1 prototype telescope has already seen first light 
ll45ll . and the Palomar Transient Factory akeady reported their 
first results 
l|46h . LSST i 
STARRS. 

These synoptic surveys will repeatedly scan the sky with re- 
visit times ("cadences") of ^ few days. The cadence timescale 
is ideally suited for following supernova light curves and de- 
tecting events near maximum brightness. Indeed, the SNLS 
have reported 289 confirmed Type la events and 117 con- 



Large-scale planned surveys include DBS 
,|4l. Sky Mapper IH], and the full-scale Pan- 



firmed core-collapse supemovae out to z '-^ 0.4 fT^. SDSS- 
II also reported 403 spectroscopically confirmed events ifTsIl 
(most of which were Type la), and 15 confirmed Type Hp 
events that are potentially capable of being used as standard- 
ized candles |44J. The Palomar Transient Factory has already 
found IztIi three events which are among the most luminous 
core-collapse events ever found, and which appear to be pul- 
sational pair-instability explosions of ultramassive stars. Fi- 
nally, Pan-STARRS 1 has reported its first confirmed super- 
nova f43\. 

Note that these surveys are unbiased in that they cover a 
large portion of the sky regions systematically and thus do not 
pre-select galaxy types or redshifts or luminosities for super- 
nova monitoring, whereas most of the past supernova surveys 
monitored pre-selected galaxies so that the results were bi- 
ased, though attempts have been made to correct for that. 

Currently, most of the design efforts for synoptic surveys 
focus on Type la supemovae, because these events are a cru- 
cial cosmological distance indicator at large redshifts. How- 
ever, the survey requirements for Type la supernova detection 
are also well-matched to collapse events, and therefore sur- 
veys that are tuned for Type la supemovae will automatically 
observe collapse events also. With their proposed properties, 
these surveys are expected to discover ~ lO'^ collapse events 
per year out to redshift z ^ 1 1221 Isol] . Due to the large sam- 
ple size, spectroscopic followup is unfeasible for most events, 
so photometric redshifts of the host galaxies (for which deep 
co-added fluxes will be available) or of the supemovae them- 
selves will be needed, just as in the case of Type la events 

iH. 

Lien & Fields 1 22] give detailed predictions for the super- 
nova harvest by synoptic surveys; here we summarize the key 
factors important for die DSNB. Within the 5-color SDSS 
ugriz bandpass system, the r and g bands provide the largest 
supernova harvest, due largely to high detector efficiency for 
these wavelengths. Moreover, distant intrinsically blue col- 
lapse events are redshifted into these bands. Detection of a 
supernova is done by differencing exposures of the same field 
of view. To determine if a transient is a supernova and to es- 
tablish its type, one must follow the supernova through the rise 
and fall of its light curve. Consequently the peak flux must be 
brighter than the minimum flux for point source detections, 
and following Ref. [22] we set a supernova limiting magni- 
tude mf/Jj^ = miim - l^'^s that is brighter by l^'^s than the 
single-visit point-source limit mum. Finally, for a given scan 
cadence timescale, a survey must trade off scan area Ailgcan 
and exposure depth rwlf^. Surveys with large scan area, such 
as Pan-STARRS and LSST, are planned to have survey depth 

The blue curve in the lower panel of Fig. [T] plots the ex- 
pected collapse event rate detected by LSST in r-band. One 
can see from the plot that in one year, LSST will have more 
than 100 supernova detections in all Az = 0.05 redshift bins 
out to redshift z - 0.9, and for z ~ 0.1 - 0.5, LSST wifl be 
able to detect more than 10* supemovae in each bin. Ref. ll22ll 
shows that Type Iln supemovae contribute to most of the de- 
tections for z^O.5 based on the luminosity functions pro- 
vided in Ref. ll24ll . Since this higher end of the detection red- 
shift range is highly affected by the small sample of Type Iln 
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in Ref. fl^, we also plot the black curve for reference to show 
a more conservative estimation that excludes Type Iln super- 
novae. One can see that the detection would reach z ^ 0.6 
in this case. The thickness of the blue and black curve repre- 
sent the statistical uncertainty (1/VN), which in most cases 
are thinner than the curve width because the uncertainty is 
very small due to the large number of supernovae. The full- 
sky fiducial supernova rate based on Ref. 1 10] is also plotted 
for comparison. The difference between the fiducial super- 
nova rate and the LSST detection rate is mainly due to survey 
depth (magnitude/flux limit), sky coverage and to a lesser ex- 
tent dust obscuration. 

A high precision measurement of the CSNR can therefore 
be done via direct counting of the enormous number of col- 
lapse events versus redshift. While a measurement of the 
CSNR shape will tests the consistency with results inferred 
from other methods, such as the star-formation history, the 
real power of synoptic surveys will be the high-statistics de- 
termination of the CSNR normalization. Note that this can 
in principle be determined by precision measurement of the 
CSNR at a single redshift bin, where the counts are largest. 
For a large survey like LSST, this should occur around z ^ 
0.3, which is set by the tradeoff of survey volume and limiting 
magnitudes [22]. In general, LSST is expected to probe the 
CSNR out to redshift z - 0.9 to 1/Vn - 10% statistical 
precision within one year of observation. 

As mentioned earlier, detections in the z ~ 0.5 — 0.9 range 
will be dominated by the most luminous core-coUapse events. 
In a study of the core-collapse luminosity function based on 
relatively sparse and inhomogeneously taken data, the rela- 
tively rare Type Iln events were found to be the most intrinsi- 
cally luminous |24]; and ultraluminous Type Iln events have 
been found |52]. Recent observations, including those by the 
synoptic Palomar Transient Factory and by ROTSE-IIIATexas 
Supernova Search, show that other core-collapse types can 
also lead to ultraluminous explosions; of these, the newly- 
discovered pair-instability outbursts are particularly intriguing 
and encouraging because this entire class of events has likely 
gone unnoticed until now ifTsL l25l - l27ll . There is clearly much 
more to be learned about about the bright end of the super- 
nova luminosity function. As more data of these ultralumi- 
nous events become available, the redshift reach of synoptic 
surveys will come into a much better focus. 



IV. IMPACT OF SYNOPTIC SURVEYS ON THE DSNB 

We are now in a position to assess the synoptic survey im- 
pact on the DSNB. Our viewpoint is to envision the situation 
several years from now, when synoptic surveys have been run- 
ning in earnest, and when the DSNB signal has been at last de- 
tected. Of course, real surveys will miss core-collapse events 
for a variety of reasons, yet following Ref. [22] we believe 
there is good reason to expect that these losses can be cali- 
brated, empirically or semi-empirically, and thus the absolute 
CSNR can be obtained out to z < 1; this should verify the al- 
ready well-determined shape of the cosmic star-formation rate 
in this regime. Furthermore, surveys will definitely measure 



the low-redshift normalization of the CSNR to high precision 
via direct counting. 

To be sure, it will be far from trivial to arrive at the un- 
derstanding we presuppose. There will be formidable astro- 
physical challenges in extracting from survey data the super- 
nova properties of interest, most importantly the event type, 
redshifts, and obscuration; less crucially for our purposes one 
would like as well the intrinsic luminosity. Ref. [22] discusses 
some reasons for optimism in the face of these challenges, 
and we also remind the reader that these issues are crucial not 
only for studies of the DSNB but also are central for other 
key topics in astrophysics and cosmology. Most notably, the 
problems of obtaining supernova type, redshift, and obscura- 
tion are at least as pressing (and in some respects more chal- 
lenging) when one uses supernovae as cosmological distance 
indicators and thus as probes of dark energy. Put differently, 
if survey supernovae are understood well enough to do dark 
energy cosmology, then we expect that the star-formation rate 
should be well-understood enough to give the DSNB source 
history out to z ^ 1, and the CSNR normalization to high 
precision. 

We now explore the impact of a CSNR determination of this 
kind. That is, we assume that one can use synoptic surveys to 
infer the absolute normalization and shape of the CSNR out 
to some redshift Zmax- In particular, Ref. Iq3r showed that all 
core-collapse types should be visible out to z^ax ^ 0.5, and 
the very bright Type Iln events should extend to Zmax ^ 1 
ll52i i53fl . Thus we will take the CSNR shape to be directly 
known from surveys to z = 1, and following Ref. l23\ we 
assume that the normalization will be very well-determined 
statistically, and so we will anticipate a measurement good to 
'57?.tot/'^tot — 5%; this error would be dominated by system- 
atic uncertainties at the most relevant redshifts. 

Referring again to Fig. [T] we compare the redshift reach of 
synoptic surveys with the redshift distribution of the DSNB 
sources. We see that the two are well matched. That is, within 
the detection energy range (~ 10 — 26 MeV positron energy), 
the neutrino sources peak within the redshift range of upcom- 
ing supernova surveys. Quantitatively, the detection rate is 
about 1.8 neutrinos/year within the detector energy range of 
10-26 MeV positron energy for neutrinos from all redshifts 
(i.e., z,„ax — 6). Of this total rate, events within redshift z = 1 
contribute 1.5 (87%) neutrinos/year, and events within red- 
shift z = 0.5 contribute 1.0 (54%) neutrinos/year. Our results 
are in good agreement with the numbers shown in Ref. [1(J 
and [54[. Therefore a large fraction of the observable neutri- 
nos come from events within z ^ 1, which is about the same 
redshift range as the upcoming supernova surveys. 

We thus see that using supernovae to directly infer the 
CSNR allows us to robustly predict a large fraction of the 
detectable neutrino events. A high precision measurement 
of the CSNR would therefore put a better constraint on the 
DSNB flux, which encodes knowledge of supernova neutrino 
physics. For example, one would then be able to distinguish 
the difference between neutrino models with different effec- 
tive temperatures, as demonstrated in Fig. |2] 

Figure ID plots the neutrino detection rates estimated based 
on models with different neutrino effective temperatures (T^ 
= 4, 6, 8 MeV, respectively) versus neutrino energy in 
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the observer's frame. The upper panel shows the current 
(^T^-tot /Tltot = 40% uncertainty in the cosmic supernova rate 
normaHzation. The bottom panel shows the future normaliza- 
tion uncertainty of <^7?.tot/7?.tot = 5% (dominated by system- 
atics), which would be achieved within one year observation 
of the upcoming supernova surveys. One can see that it is 
not easy to distinguish different neutrino models with the cur- 
rent 40% uncertainty. However, with a future 5% precision, 
it would be certainly possible to distinguish the differences 
between each models and therefore provide a way to study 
supernova neutrino physics by combining neutrino detections 
and supernova surveys. 
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FIG. 2: Upper Panel: Neutrino detection rate as a function of neu- 
trino observed energy, with different neutrino effective temperatures 
are plotted for comparison. The band thickness of the curves rep- 
resent a 57?.tot /7?.tot = 40% uncertainty in the current CSNR nor- 
malization. Bottom Panel: Same as the upper panel, but with a 5% 
normalization uncertainty instead, which is the uncertainty expected 
from upcoming supernova surveys with one year observations. 

Moreover, after several years of exposure, one might hope 
to attain statistics sufficient to measure the difference between 
the observed flux and the contributions from lower-redshift 
epochs sampled by survey supernovae. This difference en- 
codes a wealth of interesting physics and astrophysics. 

V. INVISIBLE SUPERNOVAE REVEALED 

The most dramatic possibility for a mismatch between the 
neutrino and optical supernova measures would reflect a real 
lack of optical explosions due to "invisible" supernovae. As 
mentioned in Section |ll] even in the context of conventional 
models there is a great uncertainty about whether stars with 



masses between 25 to 40 Mq explode or not. A Salpeter 
initial mass function dNi,/dm oc m~^-^^ ifssll . dictates that 
for collapse events in the 8 — IOOMq range, ~ 90% are 
core-collapse events (masses < 4OM0), which in our as- 
sumption make optically luminous explosions even for those 
that form black holes from fallback, and ^ 10% are direct- 
collapse events (> 40Mq) that are optically invisible, but 
have larger neutrino emission with greater total energy f ,y_tot 
and higher neutrino average temperature T^. A relatively con- 
servative case, which has recently been studied by Lunardini 
1I32I1 . would then assume that around 10% of collapse events 
failed to explode, hence one would expect that the neutrino 
flux from neutrino detectors would at least be ^ 10% higher 
than neutrino flux from supernova surveys. 

However, there remain large uncertainties in our qualita- 
tive understanding of massive star death, not to mention even 
larger quantitative uncertainties in neutrino and photon out- 
puts. If, as expected, the neutrino emission is larger for these 
events than for ordinary supernovae, then the signal increase 
in the detectors can be significantly lai-ger 161 (si E . 
Given these substantial uncertainties it is entirely possible that 
the invisible fraction is much higher than 10%. For example, 
one possible scenario is that supernovae that form black holes 
from fallback might actually belong to the invisible events cat- 
egory. Ref. 1 11] predicts the light curves of these fallback 
events with peak magnitudes around V — —13 to —15, which 
correspond to luminosities several orders of magnitude lower 
than ordinary core-collapse events. These authors also sug- 
gest that the total neutrino emission from the fallback events 
can be larger than normal supernovae lf56|] . Thus if we treat 
the fallback supernovae as invisible events with larger neu- 
trino emission, the invisible fraction will be higher than cur- 
rent estimates would suggests lf32ll . Therefore we will take the 
invisible fraction as an a priori free parameter, and explore 
constraints based on neutrinos and other observables. 

Fig. [3] shows several constraints on the visible supernova 
rate 7?.vis and invisible supernova rate 7?.invis ai z — 0. These 
constraints are estimated based on current data with the as- 
sumption that the shape of the CSNR is known, and we adopt 
the fiducial model described in Ref. UOl . Blue regions in the 
plot represent the allowed regions; the gray region represents 
the explicit exclusion from the non-observation of neutrinos; 
and white regions represent areas that are disallowed implic- 
itly, that is, they lie outside of current allow regions but are 
not banned directly based on current limits. 

One way to constrain 7?.tot is using the current observed 
cosmic star- formation rate. The ratio of massive star counts 
per unit mass into all stars depends only on the choice of ini- 
tial mass function; we take this ratio to be 0.007/Mq assum- 
ing the Salpeter Initial Mass Function (IMF) (53\. With the 
uncertainty ~ 20% in the cosmic star-formation rate normal- 
ization lUOtl . the upper and lower hmit of current star forma- 
tion rate ?A z — Q correspond to 7?.tot (0) = 1.25 ± 0.25 x 
10^^ yr~^ Mpc~'^. respectively, which set the darker blue re- 
gion in Fig. [3] Also, the present observed CSNR with ^ 40% 
uncertainty in its normalization is plotted as the light-blue re- 
gion in Fig. [3] which correspond to the value of 7?.tot(0) — 
1.25 ± 0.50 X 10-^yr-i Mpc^^. 

The DSNB limit in Fig. [3] shows the constraint estimated 
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from the current non-detection of the supernova neutrino 
background, which sets an upper bound of the total core- 
collapse supernova rate 7?.tot = T^vis +'7?-mvis- Ref. S points 
out that the upper limit on the neutrino flux set by Super-K in 
2003 corresponds to an upper limit of 2 events per year for a 
22.5 kton detector in the energy range of 18 - 26 MeV (see 
also Ref. for the temperature dependence of the Super-K 
limits in terms of flux instead of event rate). For the bench- 
mark Ti, = 4 MeV case, this limit allows the current 7?.tot to 
be 4.7 times bigger than current fiducial value if we assume 
all neutrino emission comes from visible events. On the other 
hand, the Super-K limit implies a current 7?.tot that is 0.64 
times smaller than our fiducial value if all neutrino emission 
comes from invisible events. Note that these two factors are 
not the same because there is more neutrino emission per in- 
visible event. 

The DSNB constraint has substantial uncertainties from 
both the visible and invisible supernova contributions. The 
neutrino emission from visible events depends on the neu- 
trino emission spectrum, i.e., temperature. To illustrate how 
this would change the DSNB limit, we also plotted the DSNB 
limit when assuming visible events have = 6 MeV instead 
of 4 MeV. The 6 MeV line intersects the TZvis axis at 1 .9 in- 
stead of 5.8 for the 4 MeV line. While the uncertainty in the 
neutrino emission from visible events would affect where the 
DSNB limit intersect with the 72.vis axis, the uncertainty in the 
neutrino emission from invisible events would change where 
the limit intersects with the 7?.invis axis. In this paper we adopt 
the highest-energy case for the neutrino emission from invis- 
ible events; however, if we choose the lowest-energy case in 
Ref. |31], then the limit would intersect with the 7?.invis axis 
at 4.6 and the whole region shown in Fig. |3] would be allowed 
by this limit and thus would give a weaker constraint. 

In addition to constraints based on current observational 
data, Kochanek et al. proposed new method of probing in- 
visible supernovae llssll . These authors suggested monitoring 
a million supergiants, in galaxies within 10 Mpc. Because 
the supergiant phase lasts ^ 10^ years, every year about one 
monitored supergiant will end its life. While some events will 
result in an ordinary optically bright supernovae, if any events 
lack optical outbursts - and are thus invisible by our defini- 
tion - they will simply disappear in sight. Considering that 
the local cosmic star-formation rate is about two times higher 
than the cosmic average, the lowest invisible event rate that 
predicts one disappearing event in the proposed five years ob- 
servation is around 7?.invis = 0.25 x 10^^ yr^^ Mpc^^. This 
line is shown as the horizontal line labeled as sensitivity to 
stellar disappearance in Fig. [3] 

Despite the preliminary nature of some of the constraints 
in Fig |3] several interesting trends already emerge. The al- 
lowed region for invisible supernovae is nonzero, but it is 
bounded and cannot be arbitrarily large. Future observa- 
tions will severely restrict the allowed region for visible su- 
pernovae. Obviously, the mere demonstration that 7?.invis is 
nonzero would immediately offer novel and unique insight 
into supernova physics. Moreover, any quantitative determi- 
nation of the absolute value of 7?.invis or the ratio 7?.invis/72.vis 
would give detailed insight into the explosion mechanism over 
the full range of core-collapse events. 
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FIG. 3: Summary of current and future constraints on the invisible 
supernova rate 7?.invis (i.e. the direct-collapse event rate in our as- 
sumption) and the visible supernova rate 7?.vis (i.e. the core-collapse 
event rate in our assumption). Blue regions are those allowed by 
current observed cosmic star-formation rate (CSFR) and CSNR and 
their uncertainties. The grey-lined region is disallowed based on the 
non-detection of the DSNB by Super-K with the assumption that 
= 4 MeV for visible events and 7.5 MeV (and also a higher total en- 
ergy) for invisible events. Another DSNB limit with = 6 MeV 
instead of 4 MeV is also plotted for comparison. The horizontal 
dashed line shows the sensitivity to stellar disappearance, which will 
directly probe the invisible supernova rate 15 SI . Note circles explored 
in Fig.|4]and stars in Fig.[5] Square marks a baseline shown in Fig.|4] 
and|3 

Also, Fig. [3] allows a larger invisible fraction than the 
/invis = 10% predicted from current theory. We marked sev- 
eral possible invisible fractions that we will discuss more in 
the figures below. The square represents a baseline, with in- 
visible fraction /invis — 0%. Circles mark possible /invis 
assuming the total CSNR is fixed to the fiducial number of 
7^tQ^ = 1.25 X 10"^yr-^ Mpc"3. The purple circle is the 
conservative case with /invis = 10%, and red circle marked 
the highest invisible fraction (/invis — 40%) one can reach 
with 7?.tot fixed. The coiTesponding changes in the DSNB de- 
tection are shown in Figure 21 where we see that when error 
in TZvis drops to 5%, it will become possible to tell the dif- 
ference between these three cases in the detectable neutrino 
energy range. The energy dependence of the fraction traces 
back to the higher energy of the neutrino flux from black hole 
forming supernovae. Therefore invisible events contribute a 
larger fraction of the neutrino flux at higher neutrino energy. 

Another set of key points in Fig. [3] are marked with stars. In 
choosing these points, we allow for the uncertainties in TZvis 
in order to explore even higher possible /invis values while 
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staying within current limits. If the visible event rate is fixed 
to the fiducial number of 7?.vis = 1-25 x 10"'' yr^^ Mpc"^, 
then the purple star marks the point with /invis = 10% adding 
to current fiducial TZ^is, and the red star marks the point with 
/invis — 17%, which is the highest /invis one can reach with 
7?.vis fixed. However, the visible event rate is quite uncertain 
and could fall substantially below our fiducial value. Includ- 
ing this uncertainty, the highest /invis that is allowed by cur- 
rent hmit is around the point marked by the orange star with 
/invis = 50%. Note that this point seems to lie just outside 
the DSNB constraint, however, one should keep in mind that 
the DSNB constraint is very sensible to theoretical assump- 
tion of the supernova neutrino emission and hence has its own 
uncertainty, as discussed earlier. 

The DSNB detections corresponding to the points marked 
by stars are shown in Figure |5] Note that the black curve with 
/invis = 0% represents the neutrino detections from the vis- 
ible events, and thus is the one that would be estimated by 
supernova surveys; the purple and red curves include differ- 
ent fractions of invisible events on top of the visible events, 
which represent those that would be detected by neutrino de- 
tectors. Therefore Fig. |5] illustrates how the differences be- 
tween DSNB from neutrino detectors and supernova surveys 
would encode information of the fraction of invisible events. 
Again, the band thickness in this figure indicates the expected 
5% uncertainty in 7?.vis, and it is clear that these three cases 
will be distinguishable. The DSNB detections for the very ex- 
treme case with /invis = 50% is plotted as the orange curve 
for comparison. 

A 50% invisible event fraction would lead to a significant 
difference between flux from neutrino detectors and super- 
nova surveys. We find that neutrinos due to invisible events 
within z ^ 1 would contribute around 75% of the event rate 
in the detectable energy range. For comparison, we expect 
the neutrinos associated with dust-obscured supernovae to be 
about ^ 20% of the signal. Thus, if the invisible event fraction 
approaches current limits, the neutrino census of supernovae 
should be able to rapidly and strongly point to the large contri- 
bution from these events. Additionally, an invisible event frac- 
tion of 50% could push the mass limit of the dkect-collapse 
events to as low as ^ 14 Mq with the Salpeter IMF. However, 
theories about supernova progenitors remain quite uncertain 
and therefore the lower mass limit implied by the invisible 
fraction is also not necessarily well-defined. Once the upcom- 
ing surveys put better constraints on the invisible fraction, one 
can hope to learn more about the mass limit of du^ect-collapse 
events. 



VI. ASTROPHYSICAL CHALLENGES AND PAYOFFS 

Our discussion until now has taken a point of view that by 
the time synoptic surveys are well under way, the loss of su- 
pernova detections from dust and survey depth can be cor- 
rected, either using the survey data themselves or from fol- 
lowup observations. In this section, we change our view- 
point from this optimistic, wide-ranging anticipation of future 
progress to a more restricted focus on the power of the survey- 
detected supernovae alone. 
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FIG. 4: One year neutrino detection as a function of neutrino ob- 
served energy. Three different fractions of the invisible events are 
plotted with 7?.tot fixed to the fiducial number. Curves with different 
colors correspond to the square/circles with the same color in Fig. [3] 
The band thickness of the curves represent 5% uncertainty expected 
from upcoming supernova surveys. 



For real surveys, some of the collapse events must be lost 
from detection mainly due to three factors: survey limiting 
magnitude, dust obscuration, and the invisible events without 
optical explosions. On the other hand, neutrino detection will 
be unaffected by any of these issues. Therefore, neutrino flux 
from neutrino detectors should exceed that estimated from su- 
pernova surveys. 

Supernova surveys thus provide a totally empirical, model- 
independent method to estimate the extreme lower limit to the 
DSNB by simply adding up the neutrino contribution from 
each supernova detected. The resulting lower bound to the 
DSNB flux is 



, survey SNe r/-. \ i 

47r N^[{l + Zi)e] 

(6) 

where each term in the sum is the flux contributed by each 
supernova observed in the survey, and the prefactor includes 
a correction for the fraction A57scan/47r of the sky covered 
by the survey. The fluxes depend on the luminosity distance 
Dl{z), which is fixed by precisely known cosmological pa- 
rameters. Notice that in this equation, only the neutrino en- 
ergy spectrum N^[{1 + Zi)e] depends on supernova and neu- 
trino physics. 

This "what you see is what you get" approach is robust but 
conservative. Namely, the result 0™''^°^ (e) will be an extreme 
lower bound for the DSNB flux. More detailed and quantita- 
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FIG. 5: Similar to Fig. [4] However, we allow larger numbers for 
7?.tot. Curves with different colors correspond to square/stars in 
Fig. [5] The band thickness of the curves represent 5% uncertainty 
expected from upcoming supernova surveys. 

live discussion can be found in Appendix lAl 

Once the DSNB is detected, the difference between the 
detected flux and the survey-based lower bound provides a 
unique measure of the events unseen by surveys. For exam- 
ple, it is conceivable that the survey predictions could exceed 
the DSNB detection (or upper limit!). This result would be 
very surprising and thus extremely tantalizing, as it would 
challenge our assumptions related to supernova physics and 
neutrino physics. In other words, this would mean that one or 
both terms in Eq.|6] the luminosity distance Dl and/or the su- 
pernova neutrino emission spectrum N^[{1 + Zi)e], might be 
wrong. But the physics behind Dl rests on well-established 
Friedmann-Robertson- Walker cosmology, and depends only 
on well-determined cosmological parameters. Thus a "DSNB 
deficit" would much more likely point to problems in the su- 
pernova emission spectrum iVy(e). Therefore, if the lower 
bound estimation (ff™^'^^{e) turn out to be higher than the ac- 
tual neutrino detections, we would be driven to rethink super- 
nova neutrinos in a way to substantially reduce the observable 
signal. 

The more likely and certainly more conventional expecta- 
tion is that when the DSNB is detected, its flux will be higher 
than the supernova survey lower bound (ff™^°^{e). In this 
case, the sign of the difference would be unsurprising, but the 
magnitude of the detected versus survey excess would still en- 
code valuable new information, such as the invisible fraction 
as discussed in the previous section. 

One might also hope for the possibility to combine survey 
supernovae and the DSNB to probe events that are optically 



visible but are lost due to dust obscuration; this could give in- 
sight into the nature and evolution of cosmic dust. To see how 
vcy Yvould change with different dust models, we examine 
with two extreme cases: (1) model with extremely low dust 
obscuration by assuming constant dust obscuration as those 
at local universe mentioned in Ref. l59ll : and (2) a model with 
very high dust obscuration by doubling the dust evolution with 
redshift compares to the model suggested in Ref. [59]. We 
find that with ra\l'^ = 23, the neutrino detection rate esti- 
mated from uncorrected supernova surveys changes by only 
^ 7% when comparing these two models. That is, dust mod- 
els (1) and (2) give 0.34 to 0.31 events per year, respectively. 
Therefore the neutrino detection rate estimated from super- 
nova survey is insensitive to the dust models and hence it will 
be difficult to use the DSNB to distinguish different dust mod- 
els with the expected survey precisions. 

VII. CONCLUSIONS 

With the next generation synoptic surveys coming online, 
a high precision measurement of the CSNR via direct count- 
ing will be achieved, and thus greatly reduce the uncertainty 
in the DSNB to a few percent. An interlocking set of strate- 
gies suggest themselves, by which one can leverage survey 
supernovae and the DSNB to probe neutrino physics as well 
as the astrophysics of cosmic supernovae. For example, the 
high-precision DSNB prediction based on supernova surveys 
would be able to distinguish supernova neutrino models with 
different neutrino temperatures. 

As we have shown, the z < 1 DSNB contribution com- 
prises most of signal at high energy > 10 MeV, and so a 
comparison of the high-energy predictions and observations 
would measure the amount of events unseen by surveys. One 
of the exciting possibilities is using the DSNB to probe the 
fraction of invisible events. With the current uncertainties, the 
observed cosmic star-formation rates and the CSNR already 
suggests possible ranges for the invisible fraction. Indeed, 
limits from present observables allows a substantial invisible 
events to up to ~ 50%, which is much higher than the fraction 
suggests by current supernova theories 10%). Once the 
upcoming synoptic surveys begin and provide high precisions 
on the CSNR and the cosmic star-formation rate, one can hope 
to reveal the fraction of invisible events. 

The current non-detection of the DSNB flux also limits the 
total supernova rate. However, this limit is sensitive to the 
theoretical assumptions of the total neutrino energy Eu^tot and 
neutrino temperature . Therefore the high precision of the 
DSNB prediction inferred from upcoming supernova surveys 
will make this limit stronger by providing knowledge of su- 
pernova neutrino physics. 

While it is unknown whether and to what degree truly invis- 
ible supernovae occur, it is certain that survey depth and dust 
obscuration will also hide supernovae from detections. To in- 
terpret the supernova data physically demands that we distin- 
guish between these factors. While the loss from survey depth 
is likely to be corrected by knowledge of supernova luminos- 
ity function, to entangle the degeneracy between dust obscu- 
ration and invisible events will be challenging. However, we 
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believe it is not impossible to discriminate the two. For ex- 
ample, there are observables across multiple wavelengths that 
can be used to estimate dust extinction. If we can constrain 
the amount of dust to a higher precision by combining all dif- 
ferent ways of measuring dust, then the dust effects can be 
modeled out ^. Hence, the only left main unknown would be 
the fraction of invisible events and we could learn this frac- 
tion by comparing the neutrino flux from neutrino detectors 
and supernova surveys. 

On the other hand, even without any extrapolations to 
the original observational data, precision measurement of the 
CSNR will be achieved by upcoming surveys, and thus will 
infer a robust lower limit of the DSNB flux by simply adding 
up the neutrino contribution from each supernova. 

We conclude by again underscoring the happy accidents 
that large-scale synoptic sky surveys will come online just at 
the time that large neutrino experiments should first discover 
the DSNB, and that the redshift reach of the two are com- 
parable. By exploiting the interconnections among the results 
from these observatories, we have a real hope of shedding new 
light into particle physics and particle astrophysics. 
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Appendix A: Surveys Set a Model-Independent Lower Bound to 
the DSNB 

As mentioned in Section |VI] a conservative and robust 
lower bound of the DSNB flux can be predicted by upcom- 
ing supernova surveys. Figure |6] shows our estimations for 
the lower bounds to the neutrino flux inferred from the core- 
collapse events detected in the r-band by a synoptic sur- 
vey. We keep A57scan fixed for simplicity, but show depen- 
dence on m\f^^ to illustrate the sensitivity to this parameter 
Planned surveys have sophisticated scan strategies using a 
variety of cadences; for reference, the largest scan areas of 
Pan-STARRS and LSST are planned to have a sensitivity of 
"^lim ~ ZS™'^^ in the bandpasses of interest. 



^ A possible cross-check here are Type la events. These are due to an older 
stellar population than core-collapse events and thus should not be pref- 
erentially obscured in their immediate locations; however, those in spiral 
galaxies will still suffer obscuration by host-galaxy disk material that hap- 
pens to lie along the line of sight. Thus Type la obscuration and reddening 
should set lower limits to the effects suffered by core-collapse events. 



FIG. 6: Upper Panel: Extreme lower bounds to the DSNB detection 
rate obtained by summing supernovae observed by surveys with dif- 
ferent limiting magnitudes, the blue curve is the limiting magnitude 
proposed by LSST and Pan-STARRS. The blue and red curves rep- 
resent a lower limit because they apply no correction for supernovae 
that are too dim or too obscured to be seen in surveys. The two black 
curves are shown for comparison: Top black curve is the DSNB flux 
from all core-collapse events in the universe out to redshift z ~ 6. 
Second top black curve is the DSNB flux from core-collapse events 
after considering dust obscuration but with infinite survey limiting 
magnitude. Results assume = 4 MeV. Middle Panel: The in- 
tegrated DSNB detection rate, i.e. the detection rate above a certain 
antineutrino energy and integrated out to = 30 MeV. The colors in- 
dicate the same features as in the top plot. Lower Panel: The fraction 
of the DSNB detection rate from the observed core-collapse events 
over those from the total collapse events. That is, a middle-panel 
red/blue curve divided by the highest black curve. Note that in this 
figure the x-axis starts at 2 MeV because no events can be detected 
below the threshold energy of 1.8 MeV. 



The upper panel shows the predicted neutrino detection 
rate from the observed core-collapse events versus neutrino 
energy. Results for the neutrino detection rate from core- 
collapse events observed with different limiting magnitude 
(from mf^^ = 23'"^s 26"^^s) are plotted. Additionally, 
the highest black curve plots the detection rate from all core- 
collapse events within the horizon (i.e., with no limiting mag- 
nitude applied) for comparison. The second highest black 
curve, also shows the detection rate for infinite survey lim- 
iting magnitude, but shows an estimate of the effect of dust 
extinction in the host galaxy. The middle panel shows the in- 
tegrated neutrino detection rate (tfj^"^^{> e) above energy e. 
In other words, this is the energy-integrated version of the up- 
per panel. The lowest panel shows the fraction of the neutrino 
detection rate from the observed supernovae over the events 
from all supernovae in the universe, that is, the correspond- 
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ing middle-panel red/blue curve divided by the highest black 
curve. 

The difference between the two black curves in Fig.|6]gives 
an indication of the neutrino contribution from dust-obscured 
supernovae. We see that an even larger effect is the loss of 
supernovae due to finite survey limiting magnitude. Note that 
when adding dust effects and limiting magnitudes, the reduc- 
tions of detection rates are more severe at low neutrino ener- 
gies. This is because observed neutrinos are redshifted, and as 
a result, a larger portion of low-energy neutrinos come from 
higher redshift where dust obscuration is more severe and su- 
pernova apparent magnitudes are dimmer because of larger 
distances. 

The observability of this energy dependence is to be un- 
derstood in the context of the energy threshold of the neu- 
trino detectors. For example, Super-K in its present form can 
discriminate from atmospheric backgrounds, and thus detect, 
cosmic neutrinos in the 18 — 26 MeV range. If Super-K is 
enhanced with gadolinium ifisll . background rejection would 
be sufficiently improved in the 10-18 MeV range to open 



this crucial window onto the DSNB. 

One sees more directly from the lower panel what portion 
of the total neutrino events detected by neutrino detector come 
from the observed core-collapse events with certain survey 
limiting magnitudes. This panel shows that ^ 18% of the 
neutrino events detected above 10 MeV are contributed by 
core-collapse events observed by surveys with a 23™^^ lim- 
iting magnitude. 

We could thus estimate the extreme lower limit to the DSNB 
to be ss 15% of the total detection events in the 10-18 
MeV range, and w 29% of the total events in the 18 - 26 
MeV range, assuming surveys with mf/J,, ~ 23™'*^. Sur- 
veys including deeper scans will see larger fractions, e.g., ap- 
proaching K, 54% of the event rate within 18-26 MeV for 
"^i'irn ~ 25™*^^. Notice that the numbers we showed above 
might be slightly lower than the percentages read directly 
from the lower panel of Fig. |6] since the numbers above are 
integrated only through the detectable energy range to reflect 
the best of what neutrino detectors would observe, while in 
Fig.|6] the numbers are integrated out to 30 MeV. 
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